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Cosmic rays are energetic charged particles from extraterrestrial sources, with the highest-energy 
events thought to come from extragalactic sources. Their arrival is infrequent, so detection requires 
instruments with large collecting areas. In this work, we report the detection of an extremely energetic 
particle recorded by the surface detector array of the Telescope Array experiment. We calculate 
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the particle’s energy as 244 + 29 (stat.) “36 (syst.) exa-electron volts (~40 joules). Its arrival direction 
points back to a void in the large-scale structure of the Universe. Possible explanations include a large 
deflection by the foreground magnetic field, an unidentified source in the local extragalactic 
neighborhood, or an incomplete knowledge of particle physics. 


lItrahigh-energy cosmic rays (UHECRs) 
are subatomic particles from extragalac- 
tic sources with energies greater than 
1 EeV (equal to 10” eV), which is about 
a million times as high as the energy 
reached by human-made particle accelerators. 
The origins of UHECRs are thought to be re- 
lated to the most energetic phenomena in the 
Universe, such as relativistic jets and outflows 
associated with black holes, gamma-ray bursts 
and relativistic flares of active galactic nuclei, 
or large-scale accretion shocks around clusters 
of galaxies (1). Alternatively, UHECRs might 
be produced by physics beyond the standard 
model of particle physics (2-4), though this 
possibility is constrained by upper bounds on 
the flux of ultrahigh-energy (UHE) photons 
(5, 6). The acceleration mechanisms of these 
particles are also unknown. Because cosmic rays 
are charged, they are deflected along their path 
to Earth by intervening galactic and extragalac- 
tic magnetic fields, so their arrival directions 
do not necessarily point to their sources. 
Although cosmic rays with energy >100 EeV 
have been observed (7), interactions with the 
cosmic microwave background radiation (CMB) 
(8) are expected to suppress the flux of UHECRs 
above 60 EeV (9, 10). This is because interac- 
tions between UHECR protons and the CMB 
can produce pions or induce photodisintegra- 
tion of heavy nuclei. The resulting break in 
the expected energy spectrum is known as the 
Greisen-Zatsepin-Kuzmin (GZK) cutoff (9, 10). 
This cutoff limits the origins of the highest- 
energy particles detected on Earth to sources 
with maximum distances of 50 to 100 Mpc, 
which have a short-enough path length for 
UHECRs to survive passage through the inter- 
galactic medium. At these distance scales, the 
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Universe is not homogeneous: Matter is con- 
centrated in a large-scale structure (LSS) com- 
posed of galaxy clusters, superclusters, filaments, 
and sheets, separated by intergalactic voids. A 
suppression of cosmic-ray flux at the highest 
energies consistent with the GZK cutoff has 
been observed (77-13). However, UHECRs with 
energies above the GZK cutoff are expected to 
be deflected less strongly by magnetic fields 
because of their higher kinetic energies, so their 
arrival directions are expected to be more closely 
correlated with their sources. 


The Telescope Array experiment 


At energies greater than 100 EeV, the flux of 
cosmic rays is less than one particle per cen- 
tury per square kilometer (72). This low flux 
can only be measured by an instrument with 
a collecting area of ~1000 km”. The energy, 
mass, and arrival direction of UHECRs can be 
reconstructed from the cascades of second- 
ary particles [an extensive air shower (EAS)] 
produced by their interaction with Earth’s 
atmosphere. Arrays of detectors, such as plas- 
tic scintillators or water-Cerenkov stations, can 
measure EASs when they reach the ground. 
The Telescope Array (TA) experiment is a 
cosmic-ray detector located in Utah, USA, at 
39.30° north, 112.91° west and 1370 m above 
sea level. It consists of a surface detector (SD) 
array with 507 stations arranged in a square 
grid. Each detector has two 3-m? layers of plastic 
scintillator that detect charged EAS particles. 
The stations are spaced by 1.2 km, giving a 
total effective area of 700 km? (14). The time- 
dependent response of the surface detectors 
is continuously monitored and calibrated by 
the detection of penetrating muons and elec- 
trons (14). The sky over the SD is viewed by 
fluorescence detectors, which directly mea- 
sure photons produced by the propagation of 
an EAS through the atmosphere, providing 
a calorimetric measurement of the shower 
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tors by determining Xmax, the atmospheric 
slant depth (in grams per square centimeter) 
at which an EAS deposits most of its energy. 
The Xmax Observable is related to the mass 
composition by a statistical analysis; the TA 
cannot determine the particle mass on an 
event-by-event basis (16). The SD measurements 
carry indirect information about the primary 
composition, which is also determined on a 
statistical basis using machine learning (1%. 

The arrival direction of a cosmic-ray particle 
is determined from the relative arrival times of 
the shower front at multiple SD stations, mea- 
sured by a time-synchronized global position- 
ing system (GPS) module mounted on each 
station. The particle density measured at a 
distance of 800 m from the EAS axis, Sgopo, is 
used as the energy indicator. Sgoo is converted 
to the primary energy of the cosmic ray as a 
function of zenith angle based on Monte Carlo 
EAS simulations (78). The SD energy scale has 
been calibrated to the calorimetric energy mea- 
sured by the fluorescence detectors with a 
factor of 1/1.27 (12). The resolution of the SD is 
1.5° in arrival direction and 15% in energy (12), 
with a systematic uncertainty of 21% (19). The 
detailed analysis procedure and event recon- 
struction are described in the supplementary 
materials (20). 


Energetic particle on 27 May 2021 


An unusually high-energy event was identi- 
fied during an arrival direction analysis (27) 
of all SD data taken between May 2008 and 
November 2021. This event triggered 23 de- 
tectors at the northwest region of the TA SD. 
The lateral density distribution (fig. S1) was 
used to determine Ssoo and hence the energy 
of the primary cosmic ray. Following our 
analysis procedure (20), we determined that 
the event on 27 May 2021 had a reconstructed 
energy of 244 + 29 (stat.) +7 (syst.) EeV in the 
detector frame. This energy is ~4 x 10’ times 
as high as the ~’7-TeV protons produced by the 
Large Hadron Collider (LHC) (22). When this 
cosmic-ray particle experienced its first collision 
with a nucleon at rest in the upper atmosphere, 
the corresponding center-of-mass energy of the 
particle collision, assuming the particle was a 
proton, was ~700 TeV. 

Figure 1A shows a map of the TA SD sig- 
nals that it recorded for the high-energy cos- 
mic ray, including the footprint of the EAS on 
the TA SD. Figure 1B shows the recorded sig- 
nal size measured at each surface detector. 
Table 1 summarizes the reconstructed prop- 
erties of the event. The waveforms recorded 
by detectors at distances above 2 km contain 
many peaks from muons induced by the ha- 
dronic interactions. With so many muon com- 
ponents, the primary particle is unlikely to 
be a photon because EASs induced by photons 
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A Surface detector array of TA 
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detector has a separate y axis. Labels indicate the detector number, total signal 


in units of the minimum ionizing particle (MIP), and the distance from the shower 
axis. Thick and thin lines (mostly overlapping) are the recorded signals in the upper and lower layers of each station. Each SD is identified by a four-digit number: 
The first two digits correspond to the column of the array in which the SD is located (numbered west to east), and the second two digits correspond to the row 

(numbered south to north). Colors correspond to those in (A). UTC, coordinated universal time. 


Relative time from earliest detector [us] 


Table 1. Reconstructed properties of the high-energy event. The reconstructed energy and Sgoo are given for the high-energy particle. The arrival 
direction is given in both the observed zenith-azimuth coordinates and the derived equatorial coordinates. The azimuth angle is defined to be anticlockwise 


from the east. The event time is expressed in UTC. 


Time (UTC) 
27 May 2021 10:35:56 


Energy (EeV) 


primarily consist of electromagnetic particles. 
We applied a neural network proton-photon 
classifier, developed for photoinduced shower 
searches using the TA SD (23, 24), to this event. 
The classifier excludes a photon as the pri- 
marty particle at the 99.986% confidence lev- 
el, instead favoring a proton as the primary 
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244 + 29 (stat.) +3 (syst.) 


Sgoo (M°) 


530 + 57 


Zenith angle 
38.6 + 0.4° 


particle. However, the classifier is unable to 
distinguish between protons and heavier 
nuclei for this event because the fluorescence 
detectors were not operating at the time (owing 
to bright moonlight). 

The core position of this event was located 
1.1 km from the northwest edge of the SD (Fig. 
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Azimuth angle 
206.8 + 0.6° 


R.A. 
259.9 + 0.6° 


Dec. 
16.1 + 0.5° 


1A). We evaluate the statistical uncertainty 
of the reconstructed energy using a detector 
simulation (72) and assuming the reconstructed 
geometry and energy parameters; we find an 
energy resolution of 29 EeV for this event. 
Assuming an energy spectrum of E “S above 
100 EeV, as previously measured using the TA 
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Fig. 2. Arrival direction 


of the high-energy event Equatorial 
compared with potential coordinates 
sources. The arrival direc- P 


tion of the 27 May 2021 
high-energy cosmic-ray 
particle (black circle) on a a 
sky map in equatorial f 
coordinates. Colored circles 
indicate calculated back- 
tracked directions 
assuming two models of 
the Milky Way regular 
magnetic field, labeled 
JF2012 (31) and PT2011 
(32). For each model, 
different symbols indicate 
the directions calculated 
for four possible primary 
species: proton (P; red), 
carbon (C; purple), silicon 


(Si; green), and iron (Fe; blue). The color bar indicates the relative flux expected 
from the inhomogeneous source-density distribution in the local LSS, smeared with a 
random Milky Way magnetic field. For comparison, nearby gamma ray-emitting 
active galactic nuclei are shown with filled diamonds and nearby starburst galaxies 
with filled stars, both with sizes that scale by the expected flux (38). The closest object 
to the proton backtracked direction in a gamma-ray source catalog (34) is the active 


SD (72), the migration effect (whereby lower 
energy showers are reconstructed with higher 
energies because of the energy resolution) is 
evaluated as -3%. We include an additional 
systematic uncertainty, owing to the unknown 
primary, of -10% in the direction of lower en- 
ergies, calculated from simulations (20). There 
was no lightning or thunderstorm activity re- 
corded in the vicinity of the TA site on 27 May 
2021 (25). 


Comparison with previous events 


Previously reported extremely high-energy cosmic- 
ray events include a 320-EeV particle in 1991 (26), 
a 213-EeV particle in 1993 (27), and a 280-EeV 
particle in 2001 (28). The 1991 event was mea- 
sured using fluorescence detectors, whereas 
the 1993 and 2001 events were both detected 
using surface detector arrays. All of these events 
were recorded by detectors in the Northern 
Hemisphere. A search in the Southern Hemi- 
sphere has not identified any events with en- 
ergy greater than 166 EeV (29), although there 
is an energy scale difference between the ex- 
periments (30). Although the event that we 
have detected was measured with a surface 
detector array, the reported energy of 244 EeV 
has been normalized to the equivalent energy 
that would have been measured with the TA 
fluorescence detector and is thus directly com- 
parable to the 1991 event. This normalization 
was performed because fluorescence detectors 
provide a direct, calorimetric measurement of 
the shower energy. The unnormalized TA SD 
reconstructed energy of 309 + 37(stat.) EeV 
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(20) is more appropriate for comparison with 
the 1993 and 2001 events. 


Possible sources of the cosmic ray 


Figure 2 shows the calculated arrival direc- 
tion of the 27 May 2021 event on a sky map in 
equatorial coordinates. The arrival direction is 
not far from the disk of the Milky Way, where 
the galactic magnetic field (GMF) is strong 
enough to substantially deflect even a parti- 
cle with an energy of 244 EeV, especially if the 
primary particle is a heavy nucleus with a 
large electric charge. The map also shows eight 
possible backtracked arrival directions, which 
we calculated (20) by assuming two GMF mod- 
els (31, 32) and four possible primary particles 
(proton, carbon nucleus, silicon nucleus, or iron 
nucleus). We used the backtracking method of 
a cosmic-ray propagation framework (33) to 
determine the arrival direction for the cosmic 
ray before it entered the Milky Way. 

We compared the arrival directions with a 
catalog of gamma-ray sources (34). We found 
that the active galaxy PKS 1717+177 is located 
within 2.5° of the calculated direction for a pro- 
ton primary. PKS 1717+177 is a flaring source 
(34); flaring sources have been proposed as 
potential cosmic-ray sources (35). However, 
its distance of ~600 Mpc (corresponding to a 
redshift of 0.137) (36) is expected to be too large 
for UHECR propagation to Earth because the 
average propagation distance at an energy of 
244: EeV is calculated to be ~30 Mpc for both pro- 
ton and iron primaries (20). We therefore dis- 
favor PKS 17171+177 as the source of this event. 
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galaxy PKS 1717+177. The dotted large circle centered around (R.A., Dec.) = (146.7° 
43.2°) indicates the previously reported TA hot spot (21). The dashed horizontal line 
indicates the limit of the TA field of view (FoV). The dotted circle centered around 
(R.A., Dec.) = (279.5% 18.0°) is the location of the Local Void (40). The galactic plane 
(G.P.) and the supergalactic plane (S.G.P.) are shown as solid and dotted curves, 
respectively. The Galactic Center (G.C.) is indicated by the cross symbol. deg., degrees. 


Figure 2 also shows the relative expected 
flux from an inhomogeneous source-density dis- 
tribution following the local LSS (37), weighted 
by the expected attenuation for a 244-EeV iron 
primary and smoothed to reflect the smearing 
resulting from turbulent magnetic fields in the 
Milky Way (20). Also shown are nearby gam- 
ma ray-emitting active galactic nuclei and star- 
burst galaxies, which have been proposed as 
possible cosmic-ray sources (38, 39). The ar- 
rival direction of this event is consistent with 
the location of the Local Void, a cavity between 
the Local Group of galaxies and nearby LSS fil- 
aments (40). There are only a small number of 
known galaxies in the void, none of which are 
expected sites of UHECR acceleration. Even 
considering the range of possible GMF deflec- 
tions and primary mass, we do not identify any 
candidate sources for this event. Only in the 
JF2012 GMF model and assuming an iron 
primary does the source direction approach a 
part of the LSS populated by galaxies. This 
backtracked direction is close to the starburst 
galaxy NGC 6946, also known as the Fireworks 
Galaxy, at a distance of 7.7 Mpc (41). However, 
NGC 6946 is not detected in gamma rays, so it 
is unlikely to be a strong source of UHECRs. 

If the energy of this event was close to the 
lower bound of its uncertainties, then the av- 
erage propagation distance is longer than we 
assumed in Fig. 2, and the deflection in the 
GMF would be larger (fig. S3). This effect would 
increase the number of possible source gal- 
axies, assuming a steady source (supplemen- 
tary text). For the alternative case of transient 
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Fig. 3. Arrival directions 
of all >100-EeV cosmic 
rays. Empty circles indicate 


Equatorial 
coordinates 


the arrival directions of P 
all cosmic rays observed by ai 
TA SD over 13.5 years of 30 : 
operation that had energies f or 

>100 EeV. The background ; 


and other symbols are 

he same as in Fig. 2. No 
lustering around the 
highest-energy event (thick 
circle) is evident. 


+ 
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sources, identifying a source is complicated 
by the time delays between electromagnetic 
radiation and charged particles because of the 
additional path lengths induced by magnetic 
deflection. We therefore cannot identify any 
potentially related transient sources. 

Nevertheless, the detection of this highly 
energetic particle allows us to estimate Do, the 
distance to the closest UHECR source (supple- 
mentary text). Assuming that the particle is an 
iron nucleus injected with an initial energy 
of Eo = 10° EeV, taking into account the en- 
ergy loss length estimated by the same prop- 
agation framework used in the backtracking 
method (42), we find Do = 10.333 Mpc. Al- 
ternatively, assuming a proton primary, we 
find Do = 27.0138 Mpc. At these energies, 
the UHECR background of distant sources is 
attenuated by the energy loss length, so only 
sources from the local Universe can contrib- 
ute. We set upper limits on the deflection by 
assuming a maximum value of the turbulent 
extragalactic magnetic field B\ms ~ 1 nG and 
a 1-Mpc characteristic length scale, finding 
<20° for iron and <1° for proton. 


Distribution of other TA events 


Figure 3 shows the arrival directions for the 
28 TA SD events with energies >100 EeV ob- 
served between May 2008 and November 2021 
using the same event selection (27). The total 
exposure is 1.6 x 10* km? sr year. No clustering 
with the highest-energy event is found. The 
244-EeV event came from a different direction 
than the TA hot spot, a 3.40 excess centered at 
right ascension (R.A.) 146.79, declination (Dec.) 
43.29, that was previously identified for events 
with energies >57 EeV (21). 

Although we expected events with energies 
above 100 EeV to be clustered, the observed 
arrival directions above 100 EeV have an iso- 
tropic distribution (Fig. 3). The lack of a near- 
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by source for the 244-EeV event could be due 
to larger magnetic deflections than predicted 
by the GMF models, caused by a heavy pri- 
mary particle or stronger magnetic fields than 
in the models. Alternatively, super-GZK UHECRs 
could indicate an incomplete understanding 
of particle physics. If there are unknown types 
of primary particles that are immune to the 
interactions with the CMB, they could retain 
their energy while traveling to Earth from more- 
distant active galaxies. We cannot distinguish 
between these possibilities with the observed 
events. 


Summary and conclusions 


We detected a particle with an energy of 244 + 
29 (stat.) +3; (syst.) EeV on 27 May 2021. The 
arrival direction of this event does not align 
with any known astronomical objects thought 
to be potential sources of UHECRs, even after 
accounting for deflection by the GMF under 
various assumptions. Comparison with other 
observed events at energies above 100 EeV 
shows an isotropic distribution with no ap- 
parent clustering. 


REFERENCES AND NOTES 


1. A.M. Hillas, Annu. Rev. Astron. Astrophys. 22, 425-444 
(1984). 

2. P. Bhattacharjee, G. Sigl, Phys. Rep. 327, 109-247 (2000). 

3. F. W. Stecker, S. T. Scully, New J. Phys. 11, 085003 
(2009). 

4. V. Berezinsky, M. Kachelriess, A. Vilenkin, Phys. Rev. Lett. 79, 
4302-4305 (1997). 

5. R. U. Abbasi et al., Astropart. Phys. 110, 8-14 (2019). 

6. P. Abreu et al., Astrophys. J. 933, 125 (2022). 

7. J. Linsley, Phys. Rev. Lett. 10, 146-148 (1963). 

8. A.A. Penzias, R. W. Wilson, Astrophys. J. 142, 419-421 
(1965). 

9. K. Greisen, Phys. Rev. Lett. 16, 748-750 (1966). 

10. G. T. Zatsepin, V. A. Kuzmin, JETP Lett. 4, 78 (1966). 

11. R. U. Abbasi et al., Phys. Rev. Lett. 100, 101101 (2008). 

12. T. Abu-Zayyad et al., Astrophys. J. Lett. 768, L1 (2013). 

13. A. Aab et al., Phys. Rev. Lett. 125, 121106 (2020). 

14. T. Abu-Zayyad et al., Nucl. Instrum. Methods Phys. Res. A 689, 
87-97 (2013). 


24 November 2023 


Relative expected flux [arbitrary units] 


15. H. Tokuno et al., Nucl. Instrum. Methods Phys. Res. A 676, 

54-65 (2012). 

16. R. U. Abbasi et al., Astrophys. J. 858, 76 (2018). 

17. R. U. Abbasi et al., Phys. Rev. D 99, 022002 (2019). 

18. D. Heck, J. Knapp, J. N. Capdevielle, G. Schatz, T. Thouw, 

“CORSIKA: A Monte Carlo code to simulate extensive air 

showers,” Forschungszentrum Karlsruhe Report, FZKA-6019 

(1998). 

19. R. U. Abbasi et al., Astropart. Phys. 80, 131-140 (2016). 

20. Materials and methods are available as supplementary 

materials. 

21. R. Abbasi et al., Astrophys. J. Lett. 790, L21 (2014). 

22. L. Evans, P. Bryant, J. Instrum. 3, SO8001 (2008). 

23. O. Kalashey, |. V. Kharuk, M. Y. Kuznetsov, G. |. Rubtsov; 
Telescope Array Collaboration, in Proceedings of 37th 
International Cosmic Ray Conference—PoS(ICRC2021) (2021), 
p. 864. 

24. |. Kharuk, O. Kalashev; Telescope Array Collaboration, in 
Proceedings of 37th International Cosmic Ray Conference—PoS 
(ICRC2021) (2021), p. 384. 

25. Lightning Exporter (Vaisala, 2023); https://lightning-exporter. 
vaisala.com. 

26. D. J. Bird et al., Astrophys. J. 441, 144 (1995). 

27. N. Hayashida et al., Phys. Rev. Lett. 73, 3491-3494 
(1994). 

28. N. Sakaki et al., in Proceedings of 27th International Cosmic 
Ray Conference (2001), pp. 333-336. 

29. P. Abreu et al., Astrophys. J. 935, 170 (2022). 

30. V. Verzi, D. Ivanov, Y. Tsunesada, Prog. Theor. Exp. Phys. 2017, 

2A103 (2017). 

31. R. Jansson, G. R. Farrar, Astrophys. J. 757, 14 (2012). 

32. M. S. Pshirkov, P. G. Tinyakov, P. P. Kronberg, 

K. J. Newton-McGee, Astrophys. J. 738, 192 (2011). 

33. R. Alves Batista et al., J. Cosmol. Astropart. Phys. 2016, 038 

(2016). 

34. S. Abdollahi et al., Astrophys. J. Suppl. Ser. 247, 33 (2020). 

35. G. R. Farrar, A. Gruzinov, Astrophys. J. 693, 329-332 (2009). 

36. D. Sowards-Emmerd, R. W. Romani, P. F. Michelson, 

S. E. Healey, P. L. Nolan, Astrophys. J. 626, 95-103 (2005). 

37. M. Y. Kuznetsov, P. G. Tinyakov, J. Cosmol. Astropart. Phys. 
2021, 065 (2021). 

38. A. Aab et al., Astrophys. J. Lett. 853, L29 (2018). 

39. R. U. Abbasi et al., Astrophys. J. Lett. 867, L27 (2018). 

40. R. B. Tully et al., Astrophys. J. 676, 184-205 (2008). 

41. J. J. Eldridge, L. Xiao, Mon. Not. R. Astron. Soc. Lett. 485, 
L58-L61 (2019). 

42. O. E. Kalashev, E. Kido, J. Exp. Theor. Phys. 120, 790-797 
(2015). 

43. T. Fujii, Data of an extremely energetic cosmic ray observed by 
a surface detector array of Telescope Array experiment, 
Zenodo (2023); https://doi.org/10.5281/zenodo.8427755. 


ACKNOWLEDGMENTS 
The Telescope Array experimental site became available 
through the cooperation of the Utah School and Institutional 


4 of 5 


RESEARCH | RESEARCH ARTICLE 


rust Lands Administration (SITLA), the US Bureau of Land 
Management (BLM), and the US Air Force. We appreciate the 
assistance of the State of Utah and the Fillmore offices of 

he BLM in crafting the Plan of Development for the site. P. A. Shea 
assisted the collaboration with valuable advice and supported the 
collaboration’s efforts. The people and the officials of Millard 
County, Utah, have been a source of steadfast and warm support 
or our work, which we greatly appreciate. We are indebted to 

he Millard County Road Department for their efforts to maintain 
and clear the roads. We gratefully acknowledge the contribution 
rom the technical staffs of our home institutions and the 
allocation of computer time from the Center for High Performance 
Computing at the University of Utah. We thank R. Cady for his 
long-standing contribution to the construction and operation of the 
detector and R. Mayta for her development of the event viewer 
ool. T. Fujii acknowledges insightful and productive discussions in 
he cosmic-ray group of Kyoto University and interdisciplinary 
communications in the Hakubi Center for Advanced Research, 
Kyoto University, and in the program for the Development of Next- 
generation Leading Scientists with Global Insight (L-INSIGHT). 
Funding: The Telescope Array experiment is supported by the 
Japan Society for the Promotion of Science (JSPS) through 
Grants-in-Aid for Priority Area 431, for Specially Promoted 
Research JP21000002, for Scientific Research (S) JP19104006, 
for Specially Promoted Research JP15H05693, for Scientific 
Research (S) JP15HO5741, for Science Research (A) JP18HO3705, 


Telescope Array Collaboration, Science 382, 903-907 (2023) 


for Young Scientists (A) JPH26707011, and for Fostering Joint 
International Research (B) JP19KK0074; by the joint research 
program of the Institute for Cosmic Ray Research (ICRR), the 
University of Tokyo; by the Pioneering Program of RIKEN for the 
Evolution of Matter in the Universe (r-EMU); by the US National 
Science Foundation awards PHY-1607727, PHY-1712517, PHY-1806797, 
PHY-2012934, and PHY-2112904; by the National Research Foundation 
of Korea (2017K1A4A3015188, 2020R1A2C1008230, and 
2020R1A2C2102800); by the Ministry of Science and Higher 
Education of the Russian Federation under the contract 075-15- 
2020-778; by IISN project no. 4.4501.18; by Belgian Science Policy 
under IUAP VII/37 (ULB); and by the Simons Foundation (00001470, 
NG). The Telescope Array was partially supported by the grants 
of the joint research program of the Institute for Space-Earth 
Environmental Research, Nagoya University; the Inter-University 
Research Program of the Institute for Cosmic Ray Research of the 
University of Tokyo; and by the foundations of Ezekiel R. Dumke and 
Edna Wattis Dumke, Willard L. Eccles, George S. Eccles, and Dolores 
Doré Eccles. The State of Utah supported the Telescope Array 
through its Economic Development Board, and the University of Utah 
through the Office of the Vice President for Research. Author 
contributions: T. Fujii identified the event, performed data analyses, 
and wrote the manuscript. R. Higuchi assisted the backtracking 
calculation, T. Sako performed the Monte Carlo simulations, 

M. Y. Kuznetsov calculated the relative expected flux, and |. Kharuk 
performed the proton-gamma ray classification. J. N. Matthews, 


24 November 2023 


P. Sokolsky, G. B. Thomson, H. Sagawa, C. H. Jui, S. Ogio, Y. Tsunesada, 
S. V. Troitsky, G. |. Rubtsov, P. G. Tinyakov, J. H. Kim, K. Fujita, and 

N. Globus discussed the results and commented on the manuscript. 
Other authors contributed to the detector construction, deployment, 
long-term data-taking and maintenance, software development, 

or review of the manuscript. All authors meet the journal's authorship 
criteria. Competing interests: There are no competing interests 

to declare. Data and materials availability: Raw data for this event 
and the analysis code that we used are archived at Zenodo (43). 
License information: Copyright © 2023 the authors, some 

rights reserved; exclusive licensee American Association for the 
Advancement of Science. No claim to original US government works. 
https://www.science.org/about/science-licenses-journal-article-reuse 


SUPPLEMENTARY MATERIALS 


science.org/doi/10.1126/science.abo5095 
Telescope Array Collaboration Authors 
Materials and Methods 

Supplementary Text 

Figs. S1 to S3 

References (44-50) 


Submitted 8 February 2022; resubmitted 6 December 2022 
Accepted 19 October 2023 
10.1126/science.abo5095 


5 of 5 


An extremely energetic cosmic ray observed by a surface detector array 


Telescope Array Collaboration*, R. U. Abbasi, M. G. Allen, R. Arimura, J. W. Belz, D. R. Bergman, S. A. Blake, B. K. 

Shin, |. J. Buckland, B. G. Cheon, T. Fujii, K. Fujisue, K. Fujita, M. Fukushima, G. D. Furlich, Z. R. Gerber, N. Globus, K. 
Hibino, R. Higuchi, K. Honda, D. Ikeda, H. Ito, A. lwasaki, S. Jeong, H. M. Jeong, C. H. Jui, K. Kadota, F. Kakimoto, O. E. 
Kalashev, K. Kasahara, K. Kawata, |. Knaruk, E. Kido, S. W. Kim, H. B. Kim, J. H. Kim, J. H. Kim, |. Komae, Y. Kubota, M. 
Y. Kuznetsov, K. H. Lee, B. K. Lubsandorzhiev, J. P. Lundquist, J. N. Matthews, S. Nagataki, T. Nakamura, A. Nakazawa, 
T. Nonaka, S. Ogio, M. Ono, H. Oshima, I. H. Park, M. Potts, S. Pshirkov, J. R. Remington, D. C. Rodriguez, C. Rott, 

G. I. Rubtsov, D. Ryu, H. Sagawa, N. Sakaki, T. Sako, N. Sakurai, H. Shin, J. D. Smith, P. Sokolsky, B. T. Stokes, T. S. 
Stroman, K. Takahashi, M. Takeda, A. Taketa, Y. Tameda, S. Thomas, G. B. Thomson, P. G. Tinyakov, I. Tkachev, T. 
Tomida, S. V. Troitsky, Y. Tsunesada, S. Udo, F. R. Urban, T. Wong, K. Yamazaki, Y. Yuma, Y. V. Zhezher, and Z. Zundel 


Science 382 (6673), . DOI: 10.1126/science.abo5095 


Editor’s summary 

Cosmic rays are charged particles from space. At low energies, they mostly originate from the Sun, whereas at 

high energies, they are expected to be emitted by nearby active galaxies. The Telescope Array Collaboration now 
reports the detection of a cosmic ray event with an energy of about 240 exa—electron volts, more than a million times 
higher than that achieved by artificial particle accelerators. Such high-energy particles should experience only small 
deflections by foreground magnetic fields, but tracing back the arrival direction shows no obvious source galaxy. The 
authors suggest that the foreground magnetic fields might be stronger than expected, or there could be unknown 
particle physics at high energies. —Keith T. Smith 
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